Zon-Cohen singularity and negative inverse temperature in a trapped particle limit 
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We study a Brownian particle on a moving periodic potential. We focus on the statistical prop- 
erties of the work done by the potential and the heat dissipated by the particle. When the period 
and the depth of the potential are both large, by using a boundary layer analysis, we calculate a 
cumulant generating function and a biased distribution function. The result allows us to understand 
a Zon-Cohen singularity for an extended fluctuation theorem from a view point of rare trajectories 
characterized by a negative inverse temperature of the biased distribution function. 

PACS numbers: 05.40.-a, 05.70.Ln, 02.50.Ey 



(N 
(N 



C 

O 

> 

m 
O 

in 
o 



X 



I. INTRODUCTION 

In 1993, the fluctuation theorem was discovered 
The theorem claims a symmetry property of the fluctu- 
ation of entropy production and provides us with a deep 
understanding of nonequilibrium physics 2- 6] . The first 
verification of the theorem in real experiments was done 
by Wang et al in 2001 Q. They used a Brownian particle 
dragged by an optical tweezer, and checked the fluctua- 
tion theorem for the work done by the tweezer. See also 
Ref. [8| for the detailed analysis of the system. In the 
stationary state, the expectation values of the work done 
by the tweezer and the heat dissipated by the particle are 
equal to each other. However, Zon and Cohen [H,[l3] pre- 
dicted that the fluctuations of the work and the heat were 
different. They pointed out that the fluctuation theorem 
for the heat was modified, while the fluctuation theorem 
for the work was valid. They called the modified relation 
an extended fluctuation theorem. See Refs. [ITI - [l8j for 
recent studies of the extended fluctuation theorem. 

In order to examine the extended fluctuation theorem, 
let us consider a cumulant generating function, 



G{h) = - log (e hQ ^ 
and a biased distribution function, 



P h (x,t) = e"* G W (s(x(t) - x)e hQW 



(1) 



(2) 



where Q(t) is the accumulated heat from t = to t — t 
and x(t) is the position of the particle at time t. The 
parameter h is called a biasing field, because one may 
understand that the right-hand side of @ is an expecta- 
tion value of 5{x{t) — x) with respect to a path probability 
given by multiplying the original path probability by a 
biasing factor e hQ(t)-tG(h) _ A hardly measurable trajec- 
tory, from which we evaluate a large value of Q(t), has a 
larger weight in with h > than in the original distri- 
bution function [© with h — 0]. This indicates that the 
biasing field is related to rare trajectories. Indeed, the 
large deviation theory connects rare trajectories with bi- 
asing fields more directly [H, [20| • 



The extended fluctuation theorem was equivalent to 
a sing ularity of the cumulant generating function ([1]) 
0, lldj . When \h\ is larger than a special value h c , G(h) 
becomes singular. In this paper, we call the singularity 
a Zon-Cohen singularity. The fact that the Zon-Cohen 
singularity emerges when the biasing field is larger than 
h c indicates that the singularity is related to rare trajec- 
tories. However, the relationship between the singularity 
and the behavior of the particle in hardly measurable 
trajectories is still unclear. Since the significance of fluc- 
tuation in nonequilibrium physics has been recognized 
recently, it is important to study the relationship using 
a systematic method and investigate whether the same 
kind of singularity occurs not only for the heat but also 
for the other quantities. 

In this paper, we consider a Brownian particle on a 
moving periodic potential. The model is the overdampcd 
case of the model studied by Lebowitz and Spohn in Ref. 
6]. By using a boundary layer analysis, we calculate 
a cumulant generating function and a biased distribu- 
tion function when the period and the depth of the po- 
tential are both large. As the result, we find that the 
biased distribution function becomes a canonical distri- 
bution function, where the inverse temperature is mod- 
ified by h. When \h\ > h c , the inverse temperature be- 
comes negative and the two limiting operations, which 
are the trapped particle limit and a limit of large ob- 
servation time, become non-interchangeable. This non- 
interchangeability corresponds to the Zon-Cohen singu- 
larity. We also check a conditional distribution function 
given Q(t)/t = q. It allows us to understand how hardly 
measurable trajectories cause the singularity. The dis- 
cussion might indicate that the same kind of singularity 
exists in the other quantities. 

The organization of the paper is the following. In Sec- 
tion nn we define a model and introduce a biased pro- 
cess. In Section IIIII and IIV1 we state main results and 
derive them. Finally, Section [V] is devoted to conclud- 
ing remarks. The Boltzmann constant fee is set to unity 
throughout the paper. 
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II. SET UP 
A. Model 

We consider a one-dimensional Brownian particle. The 
temperature of the solvent is denoted by T. The posi- 
tion of the particle is denoted by x{t) £ K. A force 
—dU(x)/dx is exerted on the particle, where U{x) is a 
periodic potential. The period of the potential is 2L. 
That is, U(x) satisfies 



U(x) = U(x + 2nL) 



(3) 



for n — ±1, ±2, • • • . We move the potential with a con- 
stant velocity v toward the negative direction of x. The 
motion of the particle is described by the Langevin equa- 
tion 



x{t) 



ld_ 

7 dy 



U{y) 



y=x(t)+vt 




(4) 



where £(t) is the Gaussian white noise that satisfies 
(£(*)) — and = 6(t — s), and 7 is a friction 

constant. In order to make the analysis easy, we intro- 
duce a new variable y(t) as the position of the particle 
measured within a reference frame that moves with the 
periodic potential. Concretely, y(t) is defined as 



y(t) = x(t) +vt- 2nL, 



(5) 



where n is an integer determined by — L < x(t) + vt — 
2nL < L. Note that y(t) is confined to [— L,L). From 
(j4]), we obtain the Langevin equation for y(t) as 



m 



7 dy 



U(y(t)) + v 




(6) 



This system can be realized in real experiments. See, for 
example, Ref. [2l|. Recently, the system has been used 
for experimental tests of some nonequilibrium relations 

We consider periodic potentials U(x) that satisfy the 
following condition 



lim 



dU(y) 



dy 



y=YL 



for < \Y\ < 1. A harmonic potential 



amid 

and a quartic potential 



0) = ht{x - 2nL) 2 



(7) 



(8) 
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FIG. 1: (Color online) A periodic potential that satisfies ©. 
We drew (8} by setting k = 2 and L = 1. 



does not satisfy ([7]). We display §8§ in Fig. [TJ 

We consider the work W(t) done by the periodic po- 
tential. Since the periodic potential exerts the force 
—dU (y) I dy\ y=x (t)+ v t on the particle and moves with the 
constant velocity —v, we find that W(t) is calculated as 



W(t) = (-v) 



d_ 

dy 



U{y) 



y—x(t)-\-vt 



(11) 



Next, we consider the heat Q(t) dissipated by the parti- 
cle. According to Sekimoto's argument [HJ], the rate of 
the heat dissipation is evaluated as 



Q(t) = xo jx-^f^t) 



(12) 



where the multiplication o represents the Stratonovich 
interpretation [2g|. We can immediately check that the 
first law of thermodynamics is satisfied. That is, 



dt 



ii 



(W(t) - Q(t)) = U(x(t 2 ) + vt 2 )- U(x(tx) + vh). 

(13) 

We also express pT|) and (IT2^) by using y(t). The result 
is 



W(t)=v—U(y(t)), 
dy 



I2T dU(y) 
7 dy 



(14) 



o^t). (15) 



U quaTt (x) = -ki{x ~ 2nL) 4 , 

where n is an integer determined by an inequality —L< 
x — 2nL < L, are the examples that satisfy (0. We 
mention that a linear potential 



^linear (x) = kl\x - 2nL\ 



(9) 

We denote by ( } p the expectation value over the noise 
(£( s ))^Lo with an initial distribution function p(yo) = 
(S(y(0) — yo))- By using the notation, we define a joint 
distribution function by 

(10) P(y ,y,t\p) = (S(y(0)-y )5(y(t)-y)) p . (16) 
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B. Biased process and cumulant generating 
functions 



We introduce a biased process. We consider a func- 
tion f(t), which depends on the trajectory of the particle 
(y( s ))*=o- We define the expectation value of f(t) in the 
biased process by 



(fit)) 



h w .hn 



-tG(h w ,h q ,t\- 



f(t)e 



W(t)h w +Q(t)h q 



(17) 



where G(h w ,h q ,t\p) is a cumulant generating function 
defined by 



1 



G(h w ,h q ,t\p) = -log(e 



t 



r W{t)h w +Q{t)h q 



(18) 



We note that the cumulant generating function corre- 
sponds to a thermodynamic free energy according to 
the thermodynamic formalism [27| . The two parame- 
ters h w and h q are called biasing fields. When we set 
h w = h q = 0, the biased expectation value (1171) returns 
to the original expectation value and the cumulant gen- 
erating function (fI5]l becomes 0. By using (|T7|) . we define 
a biased joint distribution function as 



Ph w ,h q (yo,y,t\p) = (5(y(t)-y)S(y(0)-y )) 



h w ,h a 



(19) 



In order to analyze the large t behavior of the cumu- 
lant generating function, we define two types of functions. 
The first one is 

G BcaJeii (h w , h q ) = lim G(h w ,h q ,t\p). (20) 
The second one is 

H eK (h w ,h q \p) = lim t[G(h w ,h q ,t\p) - G sca i e d(h w , h q )] . 

(|2"01) is called a scaled cumulant generating function fl!| . 
(|2"TT) is an excess quantity of the cumulant generating 
function, which was used for the calculation of an excess 
heat in Ref. [28J. We call it an excess cumulant gen- 
erating function. Since the excess cumulant generating 
function depends on the initial distribution function p, 
we explicitly indicated it in (j2"T|) . By using (f20| and (|2T|) , 
we may express G(h w ,h q ,t\p) as 

G(h w ,h q ,t\p) ~ G scalcd {h w ,h q ) + jH cx (h w ,h q \p). (22) 

We note that the difference of the left-hand side and the 
right-hand side of (|22l) is 0(e~ at ), where a is a positive 
constant. See (l6Gl). 



C. Biased distribution function and conditional 
distribution function 

Here, we show a useful relation between the biased 
distribution function and a conditional distribution func- 
tion. Let us consider a joint distribution function of y(0), 



y(t), and Q(t)/t, which is defined by 

P(yo,y,q,t\p) 

= (S(y(0) - y )S(y(t) - y)5(Q(t)/t - q)) t 



(23) 



By using this definition in the right-hand side of (fl~9"|) . we 
obtain a relation 

iogP ,h q (yo,y,t\p) + tG(o,h q ,t\p) 

( 24 ) 



log J dgP(y ,J/,?,t|p)e* 



Then, we define a function I(yo, y, q, t\p) by 



I(Vo,y,q,t\p) = --logP(y ,y,q,t\p). (25) 

For I(yo,y, q,t\p), we assume the following asymptotic 
form 

I(y , y, q, t\p) = I {q) + jh{y , y, q\ P ) + o(l/t) (26) 



when t is large. (|26l) corresponds to a large devia- 
tion property, and Io(q) is a large deviation function for 
Q(t)/t. By using the asymptotic form and a saddle point 
method, we calculate the right-hand side of (|2~3| . The 
result is 

tm^[h q q-I {q)]-h(y ,y,q*\p) + ^o{l/i), (27) 

where q* is defined as 

q* = argmax [h q q - I (q)] . (28) 
i 

From ([22]), d24j), and ([2?]>. we thus obtain 

G SC aicd(0, h q ) = max [h q q - I (q)} , (29) 

9 

and 

log Po,h q (yo , y, t \p) + i?ex(0, h q \p) 



1 



-h(yo,y,q*\p) + ^o(i/t). 



(30) 



Here, ([2Ut is a well-known relation between a large devi- 
ation function and a scaled cumulant generating function 
[H, E3|. By combining §5U§ with ([H and noticing the 
normalization condition for PQ^ q {yo,yit\p)i we arrive at 

vi ,m p {yo^y,q*Ap) . 1 ^ /+ x 
ft,h.(w),y t *b)= P(g . |t[p) +l7i°( 1 A), (3i) 

where P(q*,t\p) is a normalization constant defined by 



P(q*,t\p) ee J dy dyP(y , y, q*, t\p). (32) 

Since P(q,t\p) is the distribution function of Q(t)/t, we 
find that the biased distribution function is nothing but 
the conditional distribution function of y(0) and y(t) 
given Q(t)/t = q*. 



4 



III. RESULTS 



We denote the stationary distribution function of y by 

u,v,p 



an inverse temperature and a velocity in another system. 
The third result is about the behavior of the biased joint 

distribution function Ph xu ,h q {yo,y,t\Psi V ) when t and 



Pst V (2/)> where the superscripts U, v, and (3 indicate L are both large. That is, 



the periodic potential, the moving velocity of the poten- 
tial, and the inverse temperature of the solvent, respec- 
tively. Then, we define a canonical distribution function 

p^'Hv) by 



1 



-PU{y)+ 1 vPy 



(33) 



Z(v,p) 

where Z(v,j3) is a normalization constant determined by 



Z(v,p) = J ^ye 



(34) 



The first result is that p^ v (y) approaches Pcaa'^iv) 
when L is large: 



U,v,8 

P s t 



(y) 



rcan 



(y) 



(35) 



for y = YL (—1 < Y < 1). The definition of the symbol 
~ is the following. For functions A(Y) and B(Y), which 
depend on L, we define A(Y) ~ B(Y) for — 1 < Y < 1 as 
\im L ^ >00 (l/L)\og[A(Y)/B(Y)} = for each fixed Y. We 
use the symbol ~ throughout the paper. 

The second result is about the scaled cumulant gener- 
ating function G sca i c d(h w ,h q ). The function always be- 
comes a quadratic function in the limit L — > oo. That 
is, 

lim Gwed(fttu) h q ) = r )v 2 (h w + h q ) + Tjv 2 (h w + h q ) 2 . 

L—too 

It should be stressed that the result is always valid when- 
ever the periodic potential satisfies (0. From ([22]) and 
(f36|) . we notice that 



lim lim G(h w , h q ,t\p) 

L— >oc t— >oo 

= r )v 2 (h w + h q ) + T-yv 2 (h w + h q y 



(37) 



for any h w and hg. The system under consideration was 
analyzed in Ref [6j by Lebowitz and Spohn. They proved 
the fluctuation theorem in this system. The theorem is 
written as 

G SC aled(hw, h q ) = G sca i c d(—h w , — /3 — hq) [«-»•—« ) (38) 

where u means that the sign of v is reversed. We will 
re-derive (|38|) in Section ITVl for the sake of completeness. 
From (l38l). we obtain 



(39) 



lim lim G(h w ,h q ,t\p) 

L— ¥00 t— >OD 

= lim Mm G(-h w ,—{)- h g ,t\p)\ v ->- v , 

L—t-oo t— too 



which can also be verified by using (|37[) directly. 

Hereafter, we focus on the case in which the initial dis- 
tribution function p(y) is equal to a stationary distribu- 
tion function p^ v (y), where /?'(> 0) and v' represent 



Ph w ,h 9 (Vo,y,t\Pst V ,P ) 



(40) 



-p^(yo)p^ !A (y) + o(e- at ) 

for y = Y L, y = YL (-1 < Y , Y < 1) with 

ft = P' - hq, (41) 



V'P' +V(h q +h w ) 
P'-hq 



fa = P + h q , 



v(J3 + h q 



P + hq 



(42) 



(43) 



(44) 



From (14"TT) and (|43l) , we notice that the inverse temper- 
atures of the canonical distribution functions in (|40|) can 
become negative values. It turns out that the excess cu- 
mulant generating function has different asymptotic be- 
haviors according to whether the inverse temperature is 
negative or not. This is the fourth result. Concretely, 
when we set L to be large, the excess cumulant generat- 



ing function H cx (h w , h q \p s ^ v ' ) satisfies 

H ex (h w ,h q \p% v ''P') = 0(U(L)) 



(45) 



for pi < or pf < 0, and 

H cx (h w ,h q \p^'^') = 0(l) (46) 
for Pi > and P f > 0. From (gSJ), and (gBJ), we find 

(47) 

for Pi < or pf < 0, and 



lim lim G(h w , h q , t\p) = oo 

t— too L—>oo 



lim lim G(h w , h q , t\p) = lim lim G{h w , h q , t\p) 

t— »oo L— >oo L— >oo t— >oo 



(48) 



for Pi > and pf > 0. (gSj) shows that the two lim- 
iting operations, which are L — > oo and t — > oo, are 
interchangeable and the symmetry property of the fluc- 
tuation theorem is satisfied when ft and Pi are both pos- 
itive. However, when pi or ft is negative, the two lim- 
iting operations become non-interchangeable. If we take 
t — > oo first, the cumulant generating function satisfies 
the fluctuation theorem (|39|) . On the other hand, if we 
take L — > oo first, the cumulant generating function di- 
verges as shown in (|47p . This divergence corresponds to 
the Zon-Cohen sing ularity @,G3. 
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A. The negative inverse temperature and the 
Zon-Cohen singularity 

By substituting the explicit expression of the scaled 
cumulant generating function (|3"5j) into (|2"5)l and (121)1) , we 
obtain a relation between q* and /i g as 



2Tjv 2 



(49) 



From (|3"Tj) . (|4"9")l . and the third result stated above, we also 
obtain an expression of the joint conditional distribution 
function. That is, 



U.v'B' 



for y = r £, y = YL (-1 < Y , Y < 1) with 



1 



^iMil^! ) ~ p v£A {yo ) p v£A {y ^ _ 0{l/t) 

(50) 



, q-jir 
2Tjv 2 



where we used (fT3|) at the last line. We evaluate the ex- 
pectation value in the right-hand side by using the tra- 
jectories discussed above. The first term is approximated 
as —h q U(L)/t. The second term can be omitted by as- 
suming that the particle moves around the bottom of the 
potential during most of the time, then suddenly climbs 
up the potential just before the time t and reaches the 
top of the potential at the time t [y(t) — y]. Thus, we 
have 



G(h w ,h q ,t\p) > —±U(L). 



(56) 



This yields (|4"T)) . We note that the case that is neg- 
ative can also be discussed by following the same argu- 
ment above. The difference is that the particle goes down 
to the bottom of the potential from the top instead of 
climbing up. From these arguments, we understand how 
hardly measurable trajectories cause the singularity. 



v'0' + {q--iv 2 )/{2T 1 v) 
P- (q- 7 v 2 )/(2T^v 2 ) ' 



2Tjv 2 



Vf = V. 



(52) 



(53) 



(54) 



IV. DERIVATION 

Here, we derive the results. In the first subsection, 
we analyze the system with L fixed. Then, in the sec- 
ond subsection, we perform a boundary layer analysis 
by considering the limit L — > oo. Finally, in the third 
subsection, we derive the main results of the paper. 



From (I5TH) , (|5T1) , and (153)) , we find that the particle tends 
to climb up the potential and to reach the top of the 
potential at time t if Q(t)/t is smaller than — 'yv 2 [or to 
climb down the potential from the top at time if Q(t)/t 
is larger than ^v 2 (2j3' / (3 + 1)]. Here, we show that one 
can obtain the singularity from these rare trajectories by 
using an intuitive argument. 

Now, let us imagine that we measure trajectories and 
evaluate G(0,h q ,t\p^ v ^ ) from the trajectories. We set 
t to be sufficiently large. Then, from (l29l) . the trajecto- 

must 



U,v',B\ 



ries required for the calculation of G(0, h q , t\p, 
satisfy Q(t)/t — q*, where q* was given by (|4"9")l . Here, 
we consider the case that f3f is negative. It indicates that 
the trajectories for the calculation of G(0,h q ,t\p^ v ) 
also satisfy y(t) = L because of the negative inverse tem- 
perature. Here, by using Jensen's inequality, we obtain 

G(h w ,h q ,t\p^ v ' ) 



> 



t 



dt 



U,v> ,&' 
Pst 



-h myit ))-U(y(0))) „.„,,,, 

t Put 

h w + h a 



(55) 



t 



dt(W(t) 



A. The method of the largest eigenvalue problem 
and the Cole-Hopf transformation 



We define an operator h by 



4lU • <p 

_ _d_ 
dy 



1 d 

T r U(y) + v ) ip 



h w v ( g^ u (y) ) V 



1 (dU(y)\ 2 
l\dy) 



1 dy 
2T 



2 u(y) 



1 dy 2 



T (dU(y)Y a 



dU(y ), 
dy 



-9 



(57) 



We denote the eigenfunctions of the operator by ip n (n — 
0, 1, 2, ...) and the corresponding eigenvalues by //„ (n = 
0,1,2,...). Here, the eigenvalues are labeled such that 
Re(/it„) < Re(/i m ) for n > m, where Re(a) is the real 
part of a. We also consider the adjoint operator of h , 



6 



which is given by 

41V -f 



+ h q 



T & 



1 (dU{y)V 
J\dyj 7 dy 



2 U(y) 



7 dy 2 



(58) 



T (dU{y)\ 2 a 2T (dU{y)\ 

We denote the eigenfunctions of by 0„ (n 

0, 1, 2, ...) and the corresponding eigenvalues by v n (n 
0, 1, 2, ...). Generally, we may set 



Vn = (/■*«)* ■ 

The largest eigenvalues of and £^ h 

do not degenerate, which indicates 

= Mo- 



(59) 

are real and 
(60) 



We also note that the eigenfunctions corresponding to 
the largest eigenvalue are real. These results come from 
the Perron- Frobenius theory. See the Appendix B of Ref. 
[29l j . The orthonormal conditions for the eigenfunctions 
are 



dy(4> n (y))*i>m(y) = <W 



(61) 



(n, m — 0, 1, 2, ...), where 5 n , m is the Kronecker 6. 
Here, we define qh m ,h q (yo,y,t\p) by 



, hq (y ,y,t\p) =e tG ^ h ^P hwtK (y ,y,t\p). (62) 



As shown in Appendix^] h ^ turns out to be the time 
evolution operator of qh w ,h q (2/0, Vi t\p)- That is, 

Q^,qh w ,h q (yo,y,t\p) = ^hlji q ■ qh w ,h«(yo,y,t\p)- (63) 

We expand qh m ,h q (yo,y,t\p) by the eigenfunctions 
(V>n(2/))nLo an d solve the time evolution equation (j6"3")l 
with the initial condition qh m ,h q (yo,y,0\p) = p(yo)Hv ~ 
yo). The result is 

oo 

qh w ,h q (yo,y,t\p) = p(y )J2(Myo)T i>n(y)e' Xnt . (64) 

Then, we consider the large t behavior of 
qh w ,h q (yot y-, t\p)- The n — term becomes domi- 
nant in the right-hand side of (|64p . By combining the 
result with the definition (I62[) . we obtain 



Ph w ,h g (yo,y,t\p) 

= e {po-G(h w ,h q ,t\p))t 

p(yo)Myo)My) + o ( e -(*>-M*«i))t) 



(65) 



Furthermore, by integrating (|65p with respect to yo and 
y, and taking the logarithm of it, we also obtain 



G(h w ,h q ,t\p) 



Mo + j log coco + O ( e -(Mo-R«(Mi))^ 



where cq and cq are defined by 



c = / dycj) (y)p(y), 



Co 



— L 



dyip (y)- 



(66) 



(67) 



(68) 



(69) 
(70) 



Here, (|69|) is a well-known result [19J, |20|. There have 
been a lot of applications in which (|69p was used. See 
Ref. [6], for example. The result ([70)1 was used for 
the calculation of an excess heat in Ref. (28j . From 
(|69l) and ([70|) . we find that the expression (1651) indi- 
cates that the biased joint distribution function becomes 
p(yo)4'o(yo)' l Po(y)^ H ° x ^ hl "'' hq ^ when t is large. Essen- 
tially the same result was discussed in Ref. [301 ]. 

Next, we use the Cole-Hopf transformation in the 
largest eigenvalue problems £^ h and , and con- 

vert the largest eigenvalue problems to a non-linear eigen- 
value problem. Here, we only see the results. We define 
a non-linear operator A4h,v by 



By comparing (I22[) with (|66[) . we arrive at 

Gscalcd(h w , hq) = MCh 

H cx (h w ,h q \p) = log c c . 



dy 2 7 



1 dU \ T d 
7 dy J 7 dy 



Then, we consider a non-linear eigenvalue problem 



M h , v ■ w h . v = K h}V , 



(71) 



(72) 



where the constant Kh. v and the periodic function 
Wh, v (y) are simultaneously determined from the bound- 
ary condition Wh, v {~L) = Wh, v (L) and the normalization 
condition 



dyw htV (y) = 0. 



We introduce a potential function of w/ l;t) (y) by 



(73) 



Vh, v (y) = - dzw h . v (z) + const. (74) 
Jo 

From these preparations, we can show following relations 



h(y) = ^exp 



h q U(y) - 



v hq +h w ,v{y) 



2T 



(75) 
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My) = -s;exp 



(h q + f3) U(y) 



/3-hq-hu 



-v(y) 



2T 



G 'scaled {hw j ^g) 



2T 



(76) 
(77) 



where the coefficients (C)*C are determined from the 
normalization condition (|61|) . The derivation of these re- 
lations is shown in Appendix [Bj We note that the similar 
arguments were presented in Refs. [29l. l3l|. 



B. Boundary layer analysis with large L limit 

Here, we evaluate the asymptotic behavior of Wh tV {y) 
and K} ltV when L is large. In this calculation, we use the 
condition ([7]). 

We use a boundary layer analysis [Hf in the non-linear 
eigenvalue problem ([72]) , First, we define Wh,v(Y) by 



Wh,v{Y) = Wh, v (LY), 



(78) 



where — 1 < Y < 1. By using uih,v(Y), we rewrite the 
left-hand side of (1721) as 



2Thv 



dU(YL) 
dy 
1 dU(YL) 
1 dy 



2 7 



1 T d 



(79) 



V W h ,v + T — 



L"/ dY 



w h . 



Now, we treat L _1 as a perturbation parameter. Since 
the coefficient in front of dwh, v /dY is 0(L _1 ), we expect 
that there exists a certain region Jb in which Wft, jt ,(K) 
changes rapidly to satisfy the periodic boundary condi- 
tion and the normalization condition. That is, 



dwh, v (Y) 



dY 



> \w h ,vi Y )\ 



(80) 



for Y £ lb- When the width of the region 1^ becomes 
in the limit L — > oo, the region is called the boundary 
layer. We assume the existence of the boundary layer in 
this problem. The basic strategy of the boundary layer 
analysis is the following: (i) constructing the solutions 
inside the boundary layer (called inner solution) and out- 
side the boundary layer (called outer solution), and (ii) 
asymptotically matching those solutions so that the con- 
tinuity and the boundary conditions are satisfied. See 
Ref. [32] for more details. In this paper, we consider 
only the outer solution and obtain the leading order of 
u>h,v(Y) by using some assumptions. We state the result 
here. The derivation is shown in Appendix [Cl 

When L is large, by utilizing the condition (J7|), we 
obtain 

u>h, v (Y) 

__ (2Tjhv -1 < Y < a- 

~ |2 7 [-v(l + Th) + (l/j)dU(YL)/dy] a_ < Y < 1 

(81) 



for hv < and 



w Kv (Y) 

= J 27 
\2T-ihv 



f2 7 [-v(l + Th) + (l/^)dU(YL)/dy\ -1 < Y < 
I 2Tihv a+ < Y < 1 



(82) 



for hv > 0. The coefficients a + and a_ 
from the condition (|73p . That is, 



arc determined 



Tva- (1 + 2Th) - 7 v+ U[L) -^ (a - L) = (83) 



and 



1V a + (1 + 2Th) + 1V - U{La+) ; U{ ~ L) = 0. (84) 



The examples of (|8Tj) for the potentials J/harmo and t7 qua rt 
are displayed as green dashed lines and purple dotted 
lines in Fig. [5] We also evaluate Wh, v (Y) by numerically 
solving (fT2"j). The numerical method is the same as the 
one used in Ref. [29| . The obtained lines are displayed 
as red solid lines in Fig. [2] The figure shows that (|8Tj) 
agrees with the numerical results. 




-40 







i /-"" J 





-0.5 



0.5 



FIG. 2: (Color online) vih,v(Y) for potentials [/harmo (2/) = 
(k/2)y 2 (left) and U qu ^t(y) = (k±/4)y A (right). Quantities 
are converted to dimensionless forms by setting 7 = T — v — 
1. We fixed k = 1, fe 4 = 1 and h = -8. We set L = 18 
for ITharmo(y) and L = 3 for U qU a,rt(y)- The green dashed 
lines and purple dotted lines are obtained from (|8ip . These 
lines correspond to 14 + 36F and —16 for (7h armo (y), and 
14 + 54F 3 and —16 for U qu!iI t(y)- The red solid lines are 
numerical results. It can be seen that Q8ip agrees with the 
numerical results. 

Next, since (|72p is valid in the region that Wh, v (Y) = 
2T^hv, we obtain Kj lyV as 



(85) 



K Kv = 2T (~/v 2 h + T 1V 2 h 2 ) , 
which leads to 
G sca \cd(h w ,h q ) = jv 2 {h w + h q ) + Tjv 2 (h w + h q ) 2 . (86) 



We also check (|85l) by comparing it with numerical re- 
sults. See Fig[3J It should be stressed that (|8"6"|) is always 
valid under general periodic potentials, as long as the 
condition ([7]) is satisfied. 
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FIG. 3: (Color online) Kh. v /(2T) for potentials t/harmo(y) = 
(k/2)y 2 , C/ quart (?/) = {k±/A)y 4 , and (7i inoar (y) = ki\y\. Quan- 
tities are converted to dimensionless forms by setting 7 = 
T — v — 1. We fixed k = 1, Ua = 1 and fei = 1 and nu- 
merically evaluated Kh,v/{2T). In each figure, dashed dotted 
line (aqua blue), dotted line (purple), and dashed line (green) 
correspond to L = 6, L = 12 and L = 18 for t/harmo, L = 1, 
L — 2 and L = 3 for t/q Uar t, and L — 5, £ = 25 and L = 50 for 
Linear ■ The red solid lines in all figures denote h + h 2 , which 
is predicted by (|85[1 . As L becomes large, Kh,v approaches 
h + h 2 for L/harmo and (7 qua rt, but does not for J7u neaf . Since 
t/harmo and C/ qU art satisfy (0 but f7u aear does not, these results 
are consistent with our formulation. 



C. Derivation of the main results of the paper 

From (EU) and JSU, we find 



which indicates 



lim a± = =pl, 

L— too 



lim Wfi,v(Y) = 2T"fhv 



(87) 



for fixed Y (-1 < Y < 1). From 

Vh q +h w ,v(y) and ^-^-^-^.-^(j/) as 



we calculate 



lim V hq+hm , v (YL)/L 



= — lim — 



dzw hq+hm Az) + const. 



(89) 



= -2T-fv(h q + h w )Y + const, 
for fixed F (— 1 < Y < 1) and 
lim V-p- h hu ,-v{y)/L 

L— ¥00 



lim — 



/ dzw-p-h q -h w .-v(z) + const. 
Jo 

2Tjv(/3 + h q + h w )Y + const. 



(90) 



for fixed F(— 1 < Y < 1). By substituting these expres- 
sions into (1751) and (1761) . we thus obtain p(yo)0o(j/o) an( i 



-00 (y) as 

p(yo)^o(yo) 

~ ^7 ex P i h qU(yo) + logp(yo) + 7« C 1 ? 
for y = Y L (-1 < y < 1) and 



(91) 



Mv) ~ ^ ex Ph (/5 + ^)t/(y)+7«(/? + h q + h w )y] 

(92) 

for y = YL (—1 < y < 1), where C and C" are constants 
determined by the normalization condition (|6ip . 

Now, we derive the results in Section IIIII We recall 
that the biased distribution function becomes the orig- 
inal distribution function when we set h w = h q = 0. 
By substituting (|92|) into ([65j) . we find that the station- 
ary distribution function of this system p^ v '^ (y) satisfies 

By recalling (|86|) . we obtain (J36j) . Furthermore, be- 
cause of ([75} and ([75]) . we have 



00 (2/) = 1pQ(y)\(h w< h q ,v)^(-h w ,-l3-h q -v)- 



(93) 



Then, by utilizing ([SU]) and (|6"9")) . we obtain the fluctua- 
tion theorem for the scaled cumulant generating function 
as 



^scaled (^11; j ^g) — G sca l e d ( /l w , /? /lq)|i;_ 



(94) 

This is ([55]). 

Here, we consider the case in which the initial distri- 
bution function is equal to the stationary distribution 

function p s £ v (y). By substituting (l9Tj) and (l92j) into 
dS5j, we obtain ([t0|) . 

Finally, we derive (j45|) and ((461). By substituting (|75|) 
and (|76[) into l|70[) , and noticing the normalization con- 
dition (ffJTj) . we obtain 



H ex (h w ,h q \p) 
= log / dy exp 

log / d//exp 

log / dy 



h q U{y) + log p{y) 



(h q + P)U(y) 



Vh q+ h w ,v(y) 



2T 

V-l3-h q -h u 



-v{y) 



2T 



x exp 



By using 

where 
H 

= loe 



-/?£%) - — (Vh, +h w ,v(y) + V-[3-h q -h w -v{y)) 

(95) 

and (1901) in this expression, we obtain 
H cx {h w ,h q \p)~H + H + Z, (96) 



rfy exp [hqll (y) + logp(y) + jv (h q + h w )y], 

(97) 
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H 



log / dy exp [- (h q + /3) U(y) + jv(f3 + h q + h w )y] , 

(98) 



and 



Z = - log J dy exp [-pU(y) + jv{/3 + 2h q + 2h w )y] . 

(99) 

Z is always 0(1): 



Z = 0(1) 



(100) 



H depends on the sign of h q + f3. It becomes O(l) for 
hg + P > 0, but 0{U{L)) for h q + (3 < 0. That is, 



H 



fO(l) ^ > -/? 
\0(C/(i)) /i g <-/3. 



(101) 



Before ending the paper, we touch on a possibility that 
the singularity takes place not only for the heat but also 
for the other quantities. We showed that the singular- 
ity occurred because of the non-interchangeability of two 
types of limits. The first one was the limit of large ob- 
servation time and the second one was the trapped parti- 
cle limit. These limits become non-interchangeable since 
there exist rare trajectories in which the particle reaches 
the top of the potential (or climbs down the potential 
from the top). Aside from the heat of trapped particle 
systems, there might exist a system in which the same 
kind of singularity appears, if the system is defined by a 
limit and the limit is non-interchangeable with the limit 
of large observation time. We would like to explore such 
a system for a deeper understanding of fluctuation in 
nonequilibrium physics. 
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H depends on the initial distribution function p(y). In 
order to see it, we introduce a parameter /3 P 



j3 p = lim 



logp(i) 



L^oo U{L) ' 



(102) 
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which can take a value from to oo. It represents an 
effective temperature of the initial distribution function. 
By using /3 P , we have 

j0(U(L)) h q >p p 

H -\o(D h q < Pp . (103) 

From (|100p , (jlOip and (|103[) , we arrive at the conclusion 

(0(U(L)) h q <-f3 
H ex (h w , h q \p) = } 0(1) -j3<h g < /3 P (104) 
[0(U(L)) h q >p p . 

By setting p = p^ v (y) in it, we obtain (j45|) and (|46|) . 



V. CONCLUDING REMARKS 

In this paper, we studied the fluctuation of the work 
and the heat for a Brownian particle on a moving periodic 
potential. We considered a trapped particle limit and dis- 
cussed the Zon-Cohen singularity for an extended fluctu- 
ation theorem. As the result, we found that a conditional 
distribution function given Q(t)/t = q, where Q(t)/t was 
the rate of the heat dissipation, became a canonical dis- 
tribution function. When q was larger (or smaller) than 
a value, the inverse temperature of the canonical distri- 
bution function became negative. This indicated that 
the particle climbed up (climbed down) the potential. It 
turned out that this behavior of the particle caused the 
Zon-Cohen singularity. 



Appendix A: Derivation of (|63p 



Here, we derive (|63|) . We consider a joint distribution 
function of y(0), y(t), W(t) and Q(t), which is defined by 



P(y ,y,W,Q,t\p) 

= p(y ) (S(y(t) - y)S(W(t) - W)S(Q(t) - Q)) 



(Al) 



yo 



From the Langevin equations ©, (TT4")) . and (TT5|) . we have 
the Fokker-Planck equation for P(yo, y, W, Q, t\p) as 



dP 



- _ Av,w,Q) p 

dt - Lpp 

iv ' W ' Q) is defined by 



(A2) 



where the Fokker-Planck operator C 



-FP 



d_ 

dy 

1 



ld_ 

1 dy 

2 



U(y) + v)ip 



FP 



" { Yy u{y) ) m* 



(dU(y)y . 



T d 2 
1 dy 2 



d_ 

d~Q 



TcP_ 

7 dy 2 



T fdU{y)V d 2 
7 V dy ) dQ 



2T d 2 
7 dQdy 



dU(y ). 
dy 



(A3) 



By multiplying (|A3I) by e wh ^+Q h i ; integrating it with 
respect W and Q, and noticing the definitions ([T9"f and 
(1521 . we obtain (53J). 
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Appendix B: Derivation of (ffgj) . ((76)1 . and ((77)1 with 
the Cole-Hopf transformation 

Here, we derive ((75]), ((76]), and ((77]) from the largest 



eigenvalue problems of £^ h and \ . Similar cal- 
culations were done in Refs. (29l [3~ij] . 

First, we consider the largest eigenvalue problem of 
r (y)t 



.(»)t 



- / h m .h„ ' 



42k • = ^o^o- 



(Bl) 



By dividing this by 0o and performing some calculations, 
we obtain 



,u , , dU T 

vo = (hq + h w )v— 1 

dy 7 



d 

— (log 00 - h q U) 



(B2) 



+ 



T d 

7 dy 



d_ 

dy 



(log </) - hqU) 



Then, we introduce a potential function Vo(y) by 

V (y) = -2T (log My) -hqU(y)). (B3) 

This transformation is called the Cole-Hopf transforma- 
tion. By substituting (|B3|) into (|B2|) and combining it 
with (l60l) and (l69l). we obtain 



9Vo 
<9y 



= 2rG sca i ed (/» u „/i,) I (B4) 



Next, we consider the largest eigenvalue problem of 



(y) 

h m ,h„ » 



C hl,h a ■ V'O = MoV'o- 



(B5) 



We divide this by ipo- Then, after some calculations, we 
obtain 

^a = - \-j; + hq + h w j {~ v )-g- 



T 

1 



— ( log^ + ( /i 9 + - ] U 



d 

dy 
18U 

7 dy 



dy 



logV'o+ \U 



(B6) 



+ 



T 9 



7 <9y 

Thus, by defining 



My) 



-2T 



log^ (2/)+ ( h q + - ) U(y) 



(B7) 



substituting it into (|B6|) . and combining it with (|69[) , we 
obtain 



M. 



j3—h w —h qi —v 




2TG BC3 ,\ e <i(h w ,h q ), (B8) 



Note that the sign of the velocity in the left-hand side 
of (|B8|) is reversed. This reflects a reversed protocol of 
moving the periodic potential. 

From ((B3|, ((B4|, dH3), and (|B8|, we obtain ((7S]), fTrJ]). 
and ((77)) . We mention that the uniqueness of the solu- 
tion of the non-linear eigenvalue problem (|72p is guaran- 
teed by the Perron-Frobenius theory, because (f72"|) can 
be rewritten as the same form as (|B1[) by following the 
same calculation from (IB4I) to (IB1I). 



Appendix C: Derivation of (|81[) and (|82[) by using 
Boundary layer analysis 



Here, we derive (I5T1) and ([52"]). We consider the outer 
solution V (Y) of ((72]) . which satisfies 



(Cl) 



where w means that the left-hand side and the right- 
hand side are of the same order of magnitude. We recall 
that ((T2"]) is a quadratic equation for Wh, v (Y). Then, by 
solving this, we obtain the following identity: 



Wh,v(Y) = 7 



± 



1 dUiYL) 
v + i >- 

7 dy 

nxm , 1 dU(YL) 

v(l + 2Th) 5 

7 dy 



R(Y) , 
(C2) 



where R(Y) is defined by 

-vHTh(l + Th) + 2K h . v /j 



R(Y) 



[v(l + 2Th) - (l/~f)dU{YL)/dyY 
-2TK^L)dw Kv ldY 
[v(l + 2Th) - (l/j)dU(YL)/dyf 



(C3) 



Here, we set u>h, v (Y) = V (Y) in the right-hand side of 
(|C3j) . From (|Clj) and ((7]), we find that the second term 
of it is negligible. Then, we assume that K/ l:V , which is 
equal to 2TG sca \ c a{h Wl h q )\ hw+hq=h , is 0(1). By combin- 
ing this assumption with ((7]), we also find that the first 
term of (|C3j) is negligible. Therefore, we omit the term 
R(Y) in (|C2p . The result leads to an expression of the 
outer solution v (Y) , 



w° h JY) = 



J2T7/W 

|2 7 [-«(1 + Th) + {lh)dU(YL)/dy] . 



(C4) 



Now, in order to connect these two solutions, we use fol- 
lowing assumptions. First, we may treat the boundary 
layer of this problem as a connection between these two 
outer solutions. Second, the number of the connecting 
points should be minimized. Third, the one of the con- 
necting points must be Y = ±1 (33j. By combining these 
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assumptions with the normalization condition (|73[) . we (1821) . 
can uniquely determine the solution Wf l v (Y) as (|8ip and 
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